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In order to investigate electronic properties of recently discovered heavy fermion superconduc- 
tors CeTIns (T=Ir and Co), we employ the relativistic linear augmented-plane-wave (RLAPW) 
method to clarify the energy band structures and Fermi surfaces of those materials. The ob- 
tained energy bands mainly due to the large hybridization between Ce 4/ and In 5p states 
well reproduce the Fermi surfaces consistent with the de Haas- van Alphen experimental results. 
However, when we attempt to understand magnetism and superconductivity in CeTIns from the 
microscopic viewpoint, the energy bands obtained in the RLAPW method are too complicated to 
analyze the system by further including electron correlations. Thus, it is necessary to prepare a 
more simplified model, keeping correctly the essential characters of the energy bands obtained in 
the band-structure calculation. For the purpose, we construct a tight-binding model for CeTIns 
by including /-/ and p-p hoppings as well as f-p hybridization, which are expressed by the 
Slater-Koster integrals, determined by the direct comparison with the band-calculation result. 
Similarity and difference between Celrlns and CeCoIns are discussed based on the obtained 
tight-binding model, suggesting a significant importance of the effect of crystalline electric field 
to understand the difference in electronic properties among CeTIns. 
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§1. Introduction 

Since the establishment of the density functional the- 
ory,EP the band-structure calculation method has been 
developed as a powerful and practical tool to analyze 
electronic properties of various kinds of materials from 
the first-principle viewpoint. Especially, due to the re- 
cent rapid developments in computational technology, 
significant improvements have been achieved in the po- 
tential of the band-structure calculation.© For instance, 
nowadays it is possible to treat the complicated system 
including one hundred atoms in the unit cell. 

On the other hand, in order to include electron cor- 
relation effects beyond the Hartree-Fock approximation, 
much efforts have been made to improve the local-density 
approximation (LDA). Along this direction, there have 
been steady improvements in approximations, such as 
the generalized gradient approximation 3 ]' and LDA+U 
method.tf However, since this is essentially a trial to in- 
clude the many-body effect with high accuracy, it is dif- 
ficult to obtain rapid improvements along this direction, 
although it is a challenging and crucial problem in the 
research field of condensed matter theory In addition, 
even if such an improved approximation has been estab- 
lished in many-body physics, there still remains another 
important step to incorporate it into the band-structure 
calculation. Although this difficulty is believed to be 
overcome in future, at least at present, it becomes a hard 
task to have significant developments immediately. 

However, numbers of target materials of the band- 
structure calculation has been rapidly increased. Typ- 



ically, strongly correlated electron systems such as tran- 
sition metal oxides and /-electron compounds have at- 
tracted much attentions in recent decades. In the re- 
search field of strongly correlated electron systems, novel 
magnetism and unconventional superconductivity have 
been the central topics both from experimental and the- 
oretical sides. In order to understand those intriguing 
phenomena, it is essentially important to include corre- 
lation effects based on the concrete knowledge on elec- 
tronic properties. 

For the purpose, it is necessary to promote a couple 
of theoretical researches in parallel with different view- 
points. Namely, one research is to analyze precisely the 
energy-band structure and Fermi surfaces by using the 
appropriate band-structure calculation techniques, in or- 
der to obtain correct information about the electronic 
properties around the Fermi energy, which can be di- 
rectly compared with the result of the angle resolved pho- 
toemission experiment. Another is to construct a sim- 
plified tight-binding model to reproduce well the energy- 
band structure around the Fermi energy, in order to in- 
clude the effect of electron correlation in this simplified 
model. We believe that those two types of researches 
should be complementary to each other in order to make 
significant progress in our understandings on novel mag- 
netism and unconventional superconductivity, although 
it is quite important to develop methods and techniques 
in each research. 

As is well known, the above framework works quite 
well in 3<i-electron materials. For instance, in high- 
T c cuprates, the result of the band-structure calcula- 
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tion based on the LDA can be reproduced well by us- 
ing the tight-binding model. Further adding short-range 
Coulomb interactions to this tight-binding model, we can 
obtain the so-called Hubbard Hamiltonian, which has 
been widely accepted as a canonical model for 3d elec- 
tron systems. Based on the Hubbard Hamiltonian, it has 
been successful to understand several anomalous proper- 
ties of high- T c cuprates. Note that this success is not due 
to the speciality of 3d electron. Also in the 4d electron 
system such as Sr2Ru04, which is_an exotic material ex- 
hibiting triplet superconductivity,© three-band Hubbard 
Hamiltonian can reproduce well the Fermi surface struc- 
ture and it becomes an appropriate microscopic model to 
investigate magnetism and superconductivity in ruthen- 
atesB 

When we turn our attentions to superconductivity in 
/-electron systems, unfortunately, our understandings 
have been almost limited in the phcnomcnological level. 
One reason is that in /-electron compounds, the band- 
structure calculation method itself should be improved, 
since the relativistic effect cannot be simply neglected in 
those materials. In general, the Coulomb interaction be- 
tween nucleus and electrons in proportion to the atomic 
number becomes very strong in ions including / elec- 
trons, indicating that the velocity of electron included in 
lanthanide and actinide ions easily becomes a fraction of 
light velocity. A simple way to take into account such 
a relativistic effect is to treat the spin-orbit interaction 
as a perturbation in the non-relativistic band-structure 
calculation method, leading to a handy method to grasp 
the feature of some /-electron compounds. However, this 
is not the first-principle method to take into account the 
relativistic effect, since the magnitude of spin-orbit in- 
teraction should be adjusted by hand. This point has 
been improved by the relativistic band- structure calcula- 
tion technique, where the Dirac equations, not the non- 
relativistic Schrodinger equations, .are directly solved to 
determined the one-electron state.LF After, the improve- 
ments by Hasegawa and co-workers,BoE£P it has been 
possible to obtain the reliable electronic band-structure 
and reproduce correctly the de Haas- van Alphen (dHvA) 
experimental results for /-electron systems. 

Due to the developments in the relativistic band- 
structure calculation techniques, we seem to recover the 
situation quite similar to that of o?-electron systems, 
namely, the tight-binding model can be constructed due 
to the comparison with the relativistic band-structure 
calculations for /-electron systems. However, we should 
note that due to the large spin-orbit coupling, the mean- 
ing of one-electron state becomes complicated in /- 
electron systems, in which the total angular momentum 
j, neither spin nor orbital degree of freedom, is the only 
good quantum number to specify the one-electron state 
at each site. Thus, it is still a non-trivial problem to 
construct the tight-binding model for /-electron systems, 
even if the relativistic band-structure calculation results 
are at hand. In this paper, we attempt to show that it 
is possible to construct the tight-binding model for /- 
electron systems by using the basis to diagonalize the 
z-component of total angular momentum j, so as to re- 
produce the result of the relativistic band-structure cal- 



culation. Then, it is suggested that such a model should 
be useful for the further microscopic discussion on mag- 
netism and superconductivity in /-electron compounds. 

As a typical example for our trial to construct the 
model, in this paper we focus on recently discovered 
Ce-based sur^ci^ud.ucting compounds CeTIns (T=Rh, 
Ir, and Co)JilBEJ which are frequently referred to 
as "Ce-115" materials. A surprising point is that 
CeCoIns exhibits the superconducting transition tem- 
perature T C =2.3K, which is the highest among yet ob- 
served for heavy fermion materials at ambient pressure. 
On the other hand, Celrln5 shows T C =0.4K which is 
much less than that of CeCoIns. Note that CeRhlns is 
antiferromagnet with the Neel temperature Tn=3.8K at 
ambient pressure, while under high-pressure it becomes 
superconducting with T C =2.1K. It is not well understood 
what is the key issue to characterize the difference in the 
ground state of these Ce-115 compounds. In order to 
clarify this point, it is necessary to establish a simplified 
model to include electron correlation effects from the mi- 
croscopic viewpoint. 

Among Ce-115 compounds, the dHvA effecLhas been 
successfully observed in Celrlns and CeCoInsJi-l'Ej' both 
of which have the huge electronic specific heat coeffi- 
cient as large as several hundreds mJ/K 2 -mol. The angu- 
lar dependence of major experimental dHvA frequency 
branches is well ..explained by a quasi two-dimensional 
Fermi surfaceJiSlij) which is a clear advantage when we 
construct a model Hamiltonian, since we can restrict our- 
selves in the two-dimensional case. This is another rea- 
son why we choose the Ce-115 compound as a typical 
example to attempt to construct a tight-binding model. 

In this paper, first we employ the relativistic band- 
structure calculation method for CeTIns (T=Ir and Co) 
to investigate the electronic properties in detail. It is 
found that there exists only a small difference between 
Celrlns and CeCoIns in the energy bands around the 
Fermi level and the structure of Fermi surfaces. Namely, 
within the band-structure calculations, those two com- 
pounds cannot be distinguished, although T c 's are quite 
different. Such a difference should originate from the is- 
sue with small energy-scale which cannot be included in 
the band-structure calculations. Then, in order to inves- 
tigate such a small energy scale feature, we construct a 
tight-binding model including Ce 4/ and In bp electrons, 
and determine several parameters in the model by com- 
paring with the band-structure calculation results. We 
reanalyze the electronic properties of CeTIns to clarify 
that the effect of crystalline electric field (CEF) is a key 
issue to distinguish Celrlns and CeCoIns. 

The organization of this paper is as follows. In Sec. 2, 
we show the relativistic band-structure calculation re- 
sults both for Celrlns and CeCoIns, to discuss the ob- 
tained energy band structure and the Fermi surfaces. In 
Sec. 3, we construct the tight-binding model for CeTIn5 
by using the Slater-Koster integrals. Then, we make 
comparison between the relativistic band-structure cal- 
culation results and the dispersion obtained from our 
tight binding model. Finally in Sec. 4, we summarize 
this paper. Throughout this paper we use the energy 
units as kn—h=l. 
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Fig. 1. Crystal structure of CeTIns, called the HoCoGas-type 
tetragonal structure. 



Table I. Lattice constants and atomic positions of Celrlns and 
CeCoIns determined experimentally.tiP The labels for atoms are 
referred in Fig. 1. 





Celrln 5 


CeCoIn 5 


a 


4.666 A(=8.818 a.u.) 


4.612 A(=8.714 a.u.) 


c 


7.516 A(=14.205 a.u.) 


7.549 A(=14.264 a.u.) 


Ce 


(0, 0, 0) 


(0, 0, 0) 


Ir 


(0, 0, 1/2) 


(0, 0, 1/2) 


In(lc) 


(1/2, 1/2, 1/2) 


(1/2, 1/2, 1/2) 


In(4i) 


(0, 1/2, ±0.30528) 


(0, 1/2, ±0.305) 




(1/2, 0, ±0.30528) 


(1/2, 0, ±0.305) 



§2. Energy Band Calculations for CeTIns 

In order to calculate the electronic energy band struc- 
ture of cerium compounds, in general, relativistic ef- 
fects should be seriously taken into account .EP As is well 
known, in the hydrogen atom, v/c is just equal to the 
fine structure constant 1/137, which is negligibly small, 
where v is the velocity of electron and c is the light ve- 
locity. However, in the atom with the atomic number Z, 
v/c for the electron in the most inner orbital (K-shell) 
is given by Z/137, indicating that the velocity of the K- 
shell electron in the uranium atom with Z=92 becomes 
67% of the light velocity. This is simply understood as 
follows: Due to the strong Coulomb attraction enhanced 
by Z, electrons near the heavy nucleus must move with 
a high speed in order to keep their stationary motion. 
Thus, the relativistic effect on the inner electrons be- 
comes more significant in an atom with larger Z. 

In order to take into account major relativistic ef- 
fects such as the relativistic energy shifts, the relativistic 
screening effects, and the spin-orbit interaction in the 
energy band structure calculation, Loucks has derived 
a relativistic augmented plane wave (APW) method 
based on the Dirac one-electron wave equation.© Al- 



though his method was a natural extension of Slater's 
non-relativistic APW method to a relativistic theory, 
it included a couple of problems. For instance, the 
symmetrization of the basis functions was not taken 
into account and the method was not self-consistent. 
These points have been improved by Yamagami and 
Hasegawa.B Note that Andersen has first suggested a 
linealized technique, to include the rclativistic effect in 
the APW method.EJ 

Among several methods, in this paper we employ 
a relativistic linear augmented-plane-wave (RLAPW) 
method. The exchange and correlation potential is con- 
sidered within LDA, while the spatial shape of the one- 
electron potential is determined in the muffin-tin approx- 
imation. The self-consistent calculation is performed by 
using the lattice constants which are determined experi- 
mentally 

2.1 Basic properties of CeTIn§ 

Ce-115 materials are categolized into the HoCoGas- 
type tetragonal structure as shown in Fig. 1, which is 
characterized by the space group P4/mmm (No. 123) and 
D* h . Note that one molecule is contained per primitive 
cell. The lattice constants a and c as well as positions of 
all atoms in the unit cell are listed in Table I. 

The 4/ electrons in CeTIns (T=Ir and Co) are as- 
sumed to be itinerant. The iteration process for solv- 
ing the Dirac one-electron equation starts with the crys- 
tal charge density that is constructed by superposing 
the relativistic atomic charge densities for neutral atoms 
Ce([Xe]4/ 1 5d 1 6s 2 ), Co([Ar]3d 7 4s 2 ), Ir([Xe]4/ 14 5d 7 6s 2 ), 
and In([Kr]4(i 10 5p 1 5s 2 ), where [Xe], [Ar], and [Kr] sym- 
bolically indicate the closed electronic configuration for 
xenon, argon, and krypton, respectively. In the calcu- 
lation for the atoms, the same exchange and correla- 
tion potential are used as for the crystal. We assume 
that the Xe core state except the 5p 6 state for Ce, the 
Xe core state into the 4/ 14 for Ir, the Ar core state for 
Co, and the Kr core state for In are unchanged during 
the iteration. Namely, the frozen-core approximation is 
adopted for these core states in the calculation for the 
crystal. In the relativistic atomic calculation, the spin- 
orbit splitting in the Ce 4/, Ce 5d, Co 3d, Ir 5d, and In 
5p state are found to be 25 mRyd., 15 mRyd., 14 mRyd., 
90 mRyd., and 18 mRyd., respectively. Note here that 
mRyd. denotes milli-Rydberg and 1 Ryd.=13.6 eV. 

In each iteration step for the self-consistent calculation 
processes, a new crystal charge density is constructed us- 
ing eighteen k points, which are uniformly distributed in 
the irreducible 1/16 part of the Brillouin zone. At each k 
in the Brillouin zone, 431 plane waves are adopted under 
the condition \k + G|<4.0(27r/a) with G the reciprocal 
lattice vector, and angular momentum up to £ max =8 are 
taken into account. 

2. 2 Results {or Celrln 5 

First let us discuss the calculated results for Celrlns, 
as shown in Fig. 2, in which we depict the energy band 
structure along the symmetry axes in the Brillouin zone 
in the energy region from — 0.5Ryd. to l.ORyd. Note 
here that the three Ce bp and twenty-five In Ad bands in 
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Fig. 2. Energy band structure for Celrlns calculated by using the 
self-consistent RLAPW method. Ep indicate the position of the 
Fermi level. 



Table II. The number of valence electrons in the Ce APW sphere, 
the Ir APW sphere, and the In APW sphere partitioned into 
angular momentum. 





s 


P 


d 


/ 


Ce 


0.28 


6.16 


2.05 


1.25 


Ir 


0.54 


0.38 


7.33 


0.02 


In(lc) 


1.02 


0.66 


9.84 


0.01 


In(4i) 


3.74 


2.69 


39.31 


0.09 



the energy range between — l.ORyd. and — 0.6Ryd. are 
not shown in Fig. 2, since those bands are irrelevant to 
the present discussion. The Fermi level Ep is located at 
0.416 Ryd. and in the vicinity of Ep, there occurs a hy- 
bridization between the Ce 4/ and In 5p states. Above 
Ep near M point, the flat 4/ bands split into two groups, 
corresponding to the total angular momentum j=5/2 
(lower bands) and 7/2 (upper bands). The magnitude 
of the splitting between those groups is estimated as 0.4 
eV, which is almost equal to the spin-orbit splitting in 
the atomic 4/ state. 

The number of the valence electrons in the APW 
sphere is partitioned into the angular momentum and 
listed in Table II. There are 8.79 valence electrons out- 
side the APW sphere in the primitive cell. Each Ce APW 
sphere contains about 1.25 electrons in the / state. We 
have obtained similar results for the number^if electrons 
in the / stale per Ce APW sphere in CeSn 3 JlZP CeNi,Ef 
and CeRu2t3' , although there is difference in the num- 
ber of electrons in the d state. We have also found that 
0.1-0.2 electrons are contained in the / state of La APW 
sphere in corresponding La compounds and these elec- 
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Fig. 3. Calculated Fermi surfaces of Celrlns for (a) 13th band 
hole sheets, (b) 14th band hole sheets, (c) 15th band electron 
sheets, and (d) 16th band electron sheets. Colors indicate the 
amount of 4/ angular momentum character on each sheet of the 
Fermi surface. Red-shift indicate the increase of the admixture 
of / electrons. The center of the Brillouin zone is set at the T 
point. 



trons do not have an atomic character but a plane-wave 
character. This means that just one 4/ electron per Ce 
atom becomes itinerant in the ground state in Celrlns 
and other paramagnetic Ce compounds. 

The total density of states at Ep is evaluated as 
N(Ep) = 134.6 states/Ryd.cell. By using this value, the 
theoretical specific heat coefficient 7b an d is estimated as 
23.3 mJ/K 2 -mol. We note that the experimental elec- 
tronic specific heat coefficient 7 cxp is 750.0 mJ/K 2 -mol. 
We define the enhancement factor for the electronic 
specific heat coefficient as A=7 cxp /7band— 1, and in the 
present case, we obtain A=31.2. The disagreement be- 
tween 7band and 7exp values is ascribed to electron corre- 
lation effect and electron-phonon interactions, which are 
not fully taken into account in the present LDA band 
theory. 

Now let us discuss the Fermi surfaces of Celrln 5 . In 
Fig. 2, the lowest twelve bands are fully occupied. The 
next four bands are partially occupied, while higher 
bands are empty. Namely, the 13th, 14th, 15th, and 16th 
bands crossing the Fermi level construct the hole or elec- 
tron sheet of the Fermi surface, as shown in Fig. 3. Note 
that the Fermi surfaces from 13th, 14th, 15th, and 16th 
bands are shown in (a) , (b) , (c) , and (d) . The Fermi sur- 
face from the 13th band consists of two equivalent small 
hole sheets centered at the X points and one hole sheet 
centered at the T point. The 14th band constructs a 
large hole sheet centered at the T point, which exhibits a 
complex network consisting of big "arms" which lie along 
the edges of Brillouin zone, as observed in Fig. 3(b). The 
15th band has two kinds of sheets, as shown in Fig. 3(c). 
One is a set of two equivalent electron-like pockets, each 
of which is centered at the R point. Another sheet in the 
15th band is a large cylindrical electron sheet which is 
centered at the M point. These electron sheets are char- 
acterized by two-dimensional Fermi surfaces. Namely, 
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Fig. 4. Energy band structure calculated for CeCoIns with the 
self-consistent RLAPW method. 



they consist of the cylindrical sheet centered at the M 
point. The 16th band constructs a small electron sheet, 
which is centered at the M point, as shown in Fig. 3(d). 
These Fermi surfaces are consistent with the previous 
results obtained by the full-potential linealized APW 
method,li-S' but here we point out a possible improvement 
on our relativistic band-structure calculation in the ap- 
pearance of the Fermi surface shown in Fig. 3(d), which 
has not been obtained in the previous band-structure 
calculations. The existence of the dHvA branch corre- 
sponding to this small Fermi surface does not contradict 
the present experimental results, but such a branch has 
not been yet identified. The number of carriers contained 
in these Fermi-surface sheets are 0.039 holes/cell, 0.624 
holes/cell, 0.626 electrons/cell and 0.037 electrons/cell in 
the 13th, 14th, 15th and 16th bands, respectively. The 
total number of holes is equal to that of electrons, which 
represents that Celrlns is a compensated metal. 

Let us turn our attention to the analysis of the angular 
momentum character of the states forming various sheets 
of the Fermi surface. Here the most important quantity 
is the amount of Ce 4/ character, which is the partial 
density of states for Ce 4/ state for each point in k- 
space on the Fermi surface. This quantity is visualized 
in Fig. 3, where the admixture of the Ce 4/ states is 
increased as the red shift in color, as shown in the scale 
diagram. The broad variation of the color from blue to 
red reflects a substantial variation of the 4/ contribution 
for different groups of states which change in the range 
from about 5% to about 85%. However, this distribution 
is different from part to part on the Fermi surface. 



5 



Table III. The number of the valence electrons in the Ce APW 
sphere, the Co APW sphere, and the In APW sphere partitioned 
into angular momentum. 





s 


P 


d 


/ 


Ce 


0.29 


6.16 


2.11 


1.28 


Co 


0.50 


0.51 


7.57 


0.02 


In(lc) 


0.99 


0.60 


9.82 


0.01 


In(4i) 


3.70 


2.70 


39.23 


0.07 




Fig. 5. Calculated Fermi surfaces of CeCoIns for (a) 13th band 
hole sheets, (b)14th band hole sheets, (c)15th band electron 
sheets, and (d)16th band electron sheets. 



2.3 Results for CeCoIn 5 

In this subsection we will discuss the electronic prop- 
erties of CeCoIns. In Fig. 4, the energy band structure 
calculated for CeCoIns is shown. Some remarks to de- 
pict this figure are the same as those in Fig. 2. First 
of all, we should note that there is no qualitaive differ- 
ence in the energy-band structure between Celrlns and 
CeCoIns. Note, however, that the position of E-p is a 
little bit shifted as £ F =0.409 Ryd. 

The occupation of the bands is also the same as that 
in Celrlns. Namely, the first twelve bands are fully 
occupied. The next four bands are partially occupied, 
while higher bands are empty. The number of the va- 
lence electrons in the APW sphere is listed for each an- 
gular momentum in Table III. There are 8.63 valence 
electrons outside the APW sphere in the primitive cell 
and each Ce APW sphere contains about 1.28 electrons 
in the / state. The total density of states is calcu- 
lated at E F as N(E F ) = 154.9 states/Ryd.cell, lead- 
ing to 7band=26.8 mJ/K 2 -mol. Note that 7 oxp =300.0 
mJ/K 2 -mol for CeCoIns, smaller than that of Celrlns, 
leading to A=10.2 small compared to that of Celrlns. 

As for the Fermi surface, the situation is quite similar 
to Celrln 5 . Namely, the 13th, 14th, 15th and 16th bands 
construct the Fermi surface, as shown in Fig. 5. When 
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we compare them with those in Fig. 3, it is remarked that 
the shapes of the Fermi surfaces are essentially the same 
as those of Celrlns, although it is possible to point out 
difference of Fermi-surface structure in the orbits running 
around the bulges of the arms, which has center at the 
X point in the 14th band, and the topology of dumbbell, 
which has centered at the R point in the 15th band. 
Note also that these Fermi surfaces reproduce well the 
previous results.Q The numbers of carriers contained 
in these Fermi-surface sheets are 0.032 holes/cell, 0.627 
holes/cell, 0.629 electrons/cell, and 0.030 electrons/cell 
in the 13th, 14th, 15th, and 16th bands, respectively. 
The total number of holes is equal to that of electrons, 
indicating that CeCoIns is a compensated metal. 

The small difference between those two compounds re- 
garding the band structure around the Fermi level seems 
to be strange, if we recall the experimental facts such 
as T c and elect ronic specific heat coefficient 'Ycxp- 

For 

Celrln 5 , T C =0.2K and 7cxp =750.0 mJ/K 2 -mol, while for 
CeCoIn 5 T C =2.3K and 7cxp =300.0 mJ/K 2 -mol. Namely, 
even though the band structure is similar to each other, 
the physical properties are quite different. This point 
is difficult to understand only within the present band- 
structrue calculation results. 

§3. Tight-Binding Analysis 

In the previous section we have focussed on the 
RLAPW band calculation results, in which the overall 
framework of the band structure around the Fermi en- 
ergy E-p is determined by the broad p band and the nar- 
row / bands. We could obtain the Fermi surfaces con- 
sistent with the experimental results for both materials, 
although we have found an unexpected small difference 
in the energy bands between Celrlns and CeCoIns. 

In order to shed light on the problem, we consider the 
simplified effctive model for CeTIns by using the tight- 
binding method. For the purpose the two-dimensional 
lattice composed of Ce and In ions, are shown in Fig. 6, 
since the Fermi surfaces exhibit two-dimensionality, as 
has explained in the previous section. When we consider 
the tight-binding model on the two-dimensional lattice, 
it is enough to take 4/ orbitals of Ce and hp orbitals of 
In ions. In this sense, the model which we will construct 
is called the f-p model. 

3.1 f-p model 

As easily deduced from Fig. 6, the f-p model Hamil- 
tonian H should be written as 



H = H t + H n 



H 



fp 



H, 



CEF i 



(3.1) 



where the first and second term indicate the hopping of 
/ and p electrons, respectively, while the third term de- 
notes the f-p hybridization. The fourth term includes 
the effect of CEF. In the following we will consider each 
term. Note that in each term of H, we simply con- 
sider the nearest-neighbor hoppings of /- and p-electrons 
through a bond, since our purpose here is a construction 
of a minimal tight-binding model to discuss magnetism 
and superconductivity by further adding Coulomb inter- 
actions. 




Ce4/ 



Fig. 6. Two-dimensional lattice composed of Ce and In ions. 
Open and hatched circles indicate the position of Ce and In ions, 
respectively. 



First let us discuss the direct hopping process of /- 
electrons. Due to the spin-orbit coupling, the energy 
levels in Ce 3+ ion are split into j=5/2 sextet and j=7/2 
octet. Since the magnitude of this splitting is as large 
as 0.4 eV, it is enough to consider only the j=5/2 sex- 
tet to evaluate the effective hopping of / electrons. The 
sextet in j=5/2 are labelled by n=— 5/2, —3/2, • • •, 5/2, 
with [i the z-component of j, but due to the time rever- 
sal symmetry, those states are classified into three pairs 
characterized by up and down "pseudo" spins. 

When we define the second- quantized operator aj^ for 
the /i-state at site i, it is convenient to introduce the new 
operators /j with pseudospin a as follows: 



/ia| — fl i-5/2' 

for "a" orbitals, 



J ltd ~ 



J 15/2' 



AbT °i — 1/2' hbl fl il/2' 

for "b" orbitals, and 



(3.2) 



(3.3) 



Act — a i3/2> hcl — a i-3/2' 

for "c" orbitals. For the standard time reversal oper- 
ator K.=— io-yK, where K denotes the operator to take 
complex conjugate, we can easily show the relation 



/C/: = <J fx 
J Ira J It- 



(3.5) 



Note that this is the same definition for real spin. 
Now we express the hopping term for / electrons as 



H t = £ tf TT ,fl Ju 



(3.6) 



where tf TT , is the hopping amplitude of / electron be- 
tween t- and r'-orbitals in Ce 3+ ions connected by the 
vector a with a=x=(±l, 0) and a=y=(0,±l). The hop- 
ping amplitude can be expressed as 



and 



56 



56 



5 


-vlo 


V5 


a/10 


2 


-V2 


V5 


-V2 


1 
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VTo 


V5 \ 


VTo 
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f) 


x/5 


V2 





(3.8) 
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where (ffcr) is the Slater-Koster integral between / or- 
bitals through a bond.EP Details for the deruation of 
hopping amplitude will be discussed elsewhere E3 

Next let us consider the p-electron hopping. In this 
case, it is not necessary to consider the effect of spin-orbit 
interaction at In site. Thus, we can consider the hopping 
of p-electrons with real spin and orbitals. Furthermore, 
p z orbital does not contribute to the hopping process in 
the x-y plane. Thus, it is enough to consider the hopping 
among m=±l orbitals (m denotes the ^-component of 
angular momentum £) , since p z orbital is expressed only 
by the m=0 state among spherical harmonics for £=l. 

Note here that it is necessary to redefine spin and or- 
bital also for p electrons, which should be consistent with 
those for / electrons, since we need to consider the /- 
p mixing later. When we define the second-quantized 
operator cj mo . for the m-state with real spin a, we can 
introduce new operators P\ Ta as follows: 

(3.9) 



for "a" orbitals and 



_c iiT' p iai 



"l-U' 



P 



Ibi 



~ c i-iT' p ibi c iii' 



(3.10) 



for "b" orbitals. Again we can easily show the relation 



Kp- lTa = op- lT _ o 



(3.11) 



Note that this definition of pseudo spin is consistent with 
that of / electron. 

After some algebraic calculations, we can obtain the 
p-electron hopping term as 



a-r'er' 



(3.12) 



ISltt'o 



Where t^ TT , is the hopping amplitude of p-electron be- 
tween r- and T'-orbitals in In ions connected by the vec- 
tor a, given as 



,x = (PP&) ( 1 -1 

pTT ' 2 V "I 1 



and 



f 



y = 

prr' 



(ppa) 



l l 
l l 



(3.13) 



(3.14) 



where (ppcr) is the Slater-Koster integral between p or- 
bitals. 

Note here that p electron hopping is expressed by the 
basis of spherical harmonics, not the cubic harmonics. 
Of course, only for the purpose to show the p-electron 
hopping, cubic harmonics can provides simpler expres- 
sion, since non-zero hopping appears only between p x 
(p y ) orbitals along x- (y-) directions. However, when we 
include the f-p hybridization, it is convenient to express 
it by using the basis of spherical harmonics. 

Now we consider the f-p hybridization term. In the 
present definitions of pseudo spins for /- and p-electrons, 
it is not necessary to distinguish the real and pseudo 
spins in the f-p mixing term, if we pay due attentions 
to the Clebsch-Gordan coefficients. Due to the similar 
calculations to derive /-/ and p-p hopping term, we can 



obtain the hybridization term as 

Htp = _ E (K b ^/Lpi + b^ 

lbrr'cr 

where the hybridization is given by 



yb _ y = UP**) 



h.c), 



4\/7 




(3.15) 



(3.16) 



for b=[±l/2, ±1/2]. Note that b is the vector connect- 
ing neighbouring Ce and In ions. For b=[±l/2, =pl/2], 
For down spins, we easily obtain 



V h , f =V* 

r tt'\ y tt 



tt'I 



V, 



TT 1 ] " 



(3.17) 



By using the operators defined above, we can express the 
CEF term as 

HCEF = E B Tr>fLfir><, + A E 4^™' ^ 
Irr'a Ira 

where B TT i indicates the coefficients to express the effect 
of CEF and A is the energy level for p-orbitals measured 
from the / electron level. By consulting with the pa- 
per of HutchingsjEP it is easy to obtain for j=5/2 and 
tetragonal lattice as 



-Baa 

-Bbb 
B cc 



= 105° + 605°, 

= -8B% + 12CLB2, 

= -2B% - 18CLB2, 

= B ca = 12V5BI 



(3.19) 



where B^ are the so-called CEF parameters. 

All items included in the Hamiltonian have been es- 
tablished, but it is more convenient to express H in the 
momentum representation for the purpose to obtain the 
energy dispersion. After straightforward calculations of 
Fourier transformation, we finally obtain 

krr'o- 



(3.20) 



j = 3(/JV) 

-krr' 28 



where S TT i is the Kronecker delta and the /-electron en- 
ergy is given by 

5e k -V107 k VEe-^ 
-Vw lk 2e k -V2 7k 

v^ek -V27k e k 

(3.21) 

with e k =cosk x +cos k y and 7k=cos fc^— cos fey. The p- 
electron energy is given as 



-p _ 



(ppa) 



-7k 
e k 



w-^^v ~7k 

and the hybridization is written by 

(fpo) I Z ^ s k 

^krr'T = —Jf- -V6c k 

with Ck=cos(fc x /2) cos(k y /2) and Sk =sm (^a;/2) sin(fcy/2). 
Note that VL , , = V,* 



(3.22) 



(3.23) 
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3.2 Comparison with band calculation results 

In the previous subsection, the f-p model has been 
constructed by using seven parameters (ffo~), {ppo~), 
(fpcr), A, £?2, -B4, and Bf, which will be determined 
in the following conditions: 

(1) CEF parameters (B®, B®, Bf) are determined from 
the experimental result for magnetic susceptibility. 

(2) The positions of the top and bottom of tight-binding 
bands correspond to T point of the 25th and M point 
of 12th bands, respectively, in the RLAPW energy band 
structure. Note that those are mainly originating from 
the In bp state. 

(3) The energy difference A between / and p levels is 
determined to be positive due to the RLAPW band cal- 
culation results. 

(4) The / electron number per site is fixed as unity, 
since Ce 3+ ion contains one / electron. However, the 
total number of electrons is not integer. It is difficult to 
reproduce quantitatively the band-structure calculation 
results (see, for instance, Table II) for electron numbers 
by a two-dimensional situation. 

(5) We should reproduce the Fermi surfaces constructed 
by the 13th, 14th and 15th bands, since they include 
significant /-electron contribution. 

(6) Inequalities (/ fcr)<(ppa) and (ffcr)<(fpa) should 
hold, since p band is much wider than / band and f-p 
mixing is large from the RLAPW results. Note, how- 
ever, that the relation between (ppo~) and (fpo~) should 
be determined by the parameter fitting. 

First let us analyze the energy bands and Fermi 
surfaces of Celrln 5 based on the f-p model. As for 
CEF parameters, we can determine them as B®=— 15K, 



10000 



Bl=0. 52K, 



and -Bj=0.47K fronLthe experimental re- 



sults for magnetic susceptibility.E3 In Fig. 7(a), we plot 
(fpa) and (ppcr) as a function of A for (//cr)=4400K, 
as a typical example to explain the fitting procedure. 
Note that we determine the parameters in the two- 
dimensional tight-binding model by comparing with 
the three-dimensional band-structure calculation results. 
Thus, the present {f fo~) inevitably becomes large com- 
pared with the actual value in the three-dimensional 
system. We can see three regions regarding the num- 
ber of Fermi surfaces: One for A<1800 (region-I), two 
for 1800<A<2100 (region-II), and three for 2100<A 
(region-Ill and IV). Due to the condition (5), we need 
three Fermi surfaces. Note, however, that for A>2500 
(region-IV), {fpa) becomes smaller than (//cr), although 
we obtain three Fermi surfaces, which is in contradiction 
to condition (6). Thus, we discard the region-IV and 
choose only region-Ill, namely, 2100<A<2500. 

Note that /-electron number is fixed to unity from the 
condition (4) as shown in Fig. 7(b), while p electron num- 
ber n p is not fixed. However, as shown in Fig. 8(b), n p 
is almost 1.5, irrespective of A. If we consult with the 
band-structure calculation results, as shown in Table I, 
n p =0.66, since we consider the two-dimensional plane 
consisting of Ce and In(lc) ions. Although the band- 
structure calculation has been done in the three dimen- 
sional situation, the f-p model is constructed by assum- 
ing pure two dimensions. Thus, the discrepancy in n p be- 
tween the band-structure calculation and tight-binding 
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Fig. 7. (a) (/per) and (pprr) as a function of A for (//<r)=4400K. 
Note that the horizontal hashed line indicate the energy of (ffcr)- 
(b) / and p electron number per site as a function of A. 



model is not primary problem in the following. It may 
be possible to consider that the present f-p model is con- 
tacted with the "reservoir" of p electrons, namely, In(4i) 
sites, indicating that n v should not be determined only 
within the simple two-dimensional tight-binding model. 

After depicting the same figures as Fig. 7(a) for several 
values of (//o - ), we can determine some region in param- 
eter space to satisfy the conditions (l)-(6). In Fig. 8(a), 
we show a typical result for (//o-)=4400K and A=2300K 
with (ppcr)=6980K and (/pcr)=5685K. The Fermi level 
is located at the bottom of the Ce 4/ bands. The states 
near E-p have strong Ce 4/ character, although they in- 
clude some In bp admixture. There are three bands cross- 
ing Ep along the T-M lines. Among them two cross also 
Ep along the X-M direction, leading to two large cylin- 
drical sheets of Fermi surface. As shown by solid curves 
in Fig. 8(b), the Fermi surface of the tight-binding model 
from the first band consists of one hole sheet centered at 
the r point. The second band constructs a cylindrical 
electron sheet centered at the M point, as shown by the 
solid curves in Fig. 8(c). There exist tiny Fermi sur- 
faces around X-point, although such Fermi surfaces in 
the RLAPW results shown by dotted curves in Fig. 8(b), 
not in Fig. 8(c). This disagreement may be due to the 
simplification of the f-p model, but such tiny Fermi sur- 
faces are not important in the following analysis. The 
third band constructs a small cylindrical electron sheet 
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Fig. 8. (a)Energy band structures around Ep for Celrlns. The 
solid and dashed curves indicate the results for tight-binding 
and the RLAPW calculation, respectively. Adjusted RLAPW 
energy levels are denoted by the circled dots. Fermi-surface lines 
discussed here are (b)13th band hole sheets, (c)14th band hole 
sheets, and (d)15th band electron sheets. In (b)-(d), solid and 
dotted curves denote the Fermi-surface lines for the tight-binding 
model and the RLAPW results, respectively. 



Fig. 9. (a)Energy band structures around Ep for CeCoL-15. 
Meanings of all notations are the same as those in Fig. 8. Fermi- 
surface lines shown here are (b)13th band hole sheets, (c)14th 
band hole sheets, and (d)15th band electron sheets. 



sistent with_another experimental analysis for magnetic 
excitations EB 



centered at the M point, as shown in Fig. 8(d). In to- 
tal, the overall feature of energy bands around the Fermi 
level are well reproduced by the simplified f-p model. 

Now we turn our attention to CeCoIns. CEF pa- 
rameters of this compound have been determined exper- 
imentally as Bg=-6.8K, B2=0.05K, and 5|=2.5K.@> 
By using these CEF parameters, we have attempted 
to reproduce the band-structure calculations, and fi- 
nally we can obtain Fig. 9 for A=2300K, (ppcr)=5730K, 
(/pcr)=5360K, and (//<t)=4400K. Again we see that ma- 
jor features of the energy bands of CeCoIn5 are basically 
similar to those of Celrlns. The behavior of the hy- 
bridization of energy bands and the relative position of 
the 4/-bands to the other valence bands are not changed 
when we compare with those of Celrlns. In fact, the ob- 
tained Fermi surfaces are essentially the same as those 
of Celrlns , as easily understood from the comparison be- 
tween Figs. 8 and 9. 

Here let us consider the meaning of similarity between 
those two compounds. In the f-p model we have seven 
parameters, which are classified into two groups, i.e., 
high- and low-energy groups. The high-energy group 
includes (ffcr), (ppcr), (fpa), and A, while low-energy 
group consists of CEF parameters, B®, and B\. Note 
here that the typical energy scale is 1000K for the for- 
mer, while at most 10K for the latter group. The energy 
band structure with the RLAPW method around the 
Fermi energy E-p is mainly determined by high-energy 
group parameters. On the other hand, the band struc- 
ture is not sensitive to CEF parameters, since those are 
much smaller than the typical energy scale considered in 
the band-structure calculations. In fact, the resolution 
of the energy band structure with the RLAPW method 
is less than the energy scale of CEF parameters. Thus, 
the smallness of difference between Celrlns and CeCoIn 5 
suggests that CEF parameters are important to charac- 
terize the difference in Ce-115 materials. This is con- 



3.3 Role of CEF parameters in Ce-115 compounds 

Since we have confirmed that the overall band- 
structure is well reproduced by the tight-binding f-p 
model, let us now discuss the role played by the CEF 
parameters in Ce-115 materials based on the f-p model 
from the qualitative viewpoint. Note again that the 
effect of CEF parameters cannot be discussed in the 
RLAPW results, since it is included only partially in the 
band-structure calculations. In order to see the effect of 
CEF parameters on the band-structure, we consider two 
types of modification in CEF parameters for Celrlns: 
In Fig. 10(a), B% is just increased as 105° by keeping 
others invariant, while only B\ is increased as 10B® in 
Fig. 10(b). As for Bf, we show no figure, since the band- 
structure is not changed significantly, even if we change 
it as 100£?|. Namely, B\ is the most insensitive quan- 
tity among three CEF parameters, and thus, we do not 
discuss the effect of B\ furthermore. 

In Figs. 10(a) and (b), it is noted that there is clear 
difference in the tendency of the change in the band- 
structure when we increase B® or B^. Here note that 
the solid curves indicates the results for increasing CEF 
parameter. When we increase B®, the energy bands are 
totally shifted upward, as observed in Fig. 10(a). On 
the other hand, for the increase of B®, the distance be- 
tween two bands forming the Fermi surfaces tend to be 
narrow, as shown in Fig. 10(b). Here we focus on two 
bands forming the main Fermi surfaces, second and third 
solid curves in Fig. 8 from the bottom, since those Fermi 
surfaces with large volume in the Brillouin zone includes 
significant amount of / electrons. Note that this nar- 
rowing effect due to the increase of B® occurs only in 
bands around the Fermi level, although the shift due to 
the change of B% seems to occur in all bands. 

It is considered that the physical quantities are sen- 
sitive to -B4, since those are determined by the band- 
structure around the Fermi level. For instance, due 
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Fig. 10. Energy band structures of the f-p model around E-p for 
Celrlns when we increase (a) only BJJ and (b) only B^, by keep- 
ing other parameters invatiant. The Fermi level for the band- 
structure for lO-Bj ' s fixed as that for B®, and /-electron number 
is unity from the condition (4). 



binding model for quasi-particles, since we compare it 
with the band-structure calculation results, even though 
correlation effects have not been satisfactrily included in 
the band-structure calculations. In any case, in order to 
obtain definite conclusion on this issue, further quanti- 
tative works to include electron correlations are needed 
in future. 

Finally, let us try to consider a possible scenario to 
understand the difference in T c between CeCoIn5 and 
Celrlns. It may be risky to conclude something regard- 
ing superconductivity without considering seriously the 
effect of electron correlations, since superconductivity in 
Ce-115 materials has been believed to.be originating from 
antiferromgantic spin fluctuations.^ However, it is an 
interesting trial to discuss the effect of B® on T c to clarify 
the importance of CEF parameters in Ce-115 materials. 
As discussed above, the increase in induces the effec- 
tive narrowing of the band- width around the Fermi level. 
Thus, it is deduced that electron correlation becomes ef- 
fectively larger as B® is increased. In fact, for CeRhlns 
which is an antiferr.Gunagnet at ambient pressure, B® is 
estimated as 0.55KJH3) which is slightly larger than that 
of Celrlns. This is consistent with the present conclu- 
sion. By following our discussion, it may be possible 
to conclude that the magnitude of antiferromgantic spin 
fluctuations should be larger in the order of CeCoIns, 
Celrlns, and CeRhIn5, in agreement with the previous 
discussion on the effect of level-splitting to control anti- 
ferromgantic spin fluctuations Ef Of course, the present 
discussion on supercodnuctivity is just in a qualitative 
level, but we believe that the effect of CEF should be 
one of key issues to distinguish the differece in Ce-115 
materials, in combination with the large renormalization 
effect due to electron correlations. 



to the narrowing effect of effective band-width around 
the Fermi level, N(Ep) (the density of states at the 
Fermi level) should be increased with increasing B®. 
Here we recall the experimetal results on CEF param- 
eters such as B^0.52K for Celrln 5 and B^O.ObK for 
CeCoIns. Due to the present discussion, we can con- 
clude that N(Ep) of Celrln 5 becomes larger than that of 
CeCoIns, if we correctly include the effect of CEF in the 
band-structure calculations. In fact, we have observed 
that 7cxp =750.0 mJ/K 2 -mol and 300.0 mJ/K 2 -mol for 
Celrlns and CeCoIns, respectively. We emphasize that 
the difference in 7 cxp can be understood, at least qual- 
itatively, by taking into account the effect of CEF pa- 
rameters into the f-p model. It is a clear advantage to 
consider the elecrtonic properties based on the present 
tight-binding model, since the effect of CEF parameters 
can be included in an explicit manner. 

Here it should be also noted that at present, electron 
correlations have not been included explicitly. In general, 
in actual heavy-fermion materials, quantities with small 
energy scale such as CEF parameters can be significant, 
since the energy scale of quasi-particles is much reduced 
due to the large renormalization effect. In this sense, 
it may be considered that our present analysis on the 
CEF parameters has been done effectively on the tight - 



§4. Summary 

In this paper, we have applied the RLAPW method to 
the self-consistent calculation of the electronic structure 
for Celrlns and CeCoIns on the basis of the itinerant 
4/ electron picture. We have found that a hybridization 
between the Ce 4/ state and In 5p state occurs in the 
vicinity of E-p. The obtained main Fermi surfaces are 
composed of two hole sheets and one electron sheet, all 
of which are constructed from the band having the Ce 
4/ state and the In bp state. 

In order to understand the difference in electronic 
properties between Celrlns and CeCoIns even though 
the energy-band structure around the Fermi level are 
quite similar to each other, we have constructed the /- 
p model based on the knowledge of the RLAPW band 
calculation. It has been found that the size and position 
of the main Fermi surfaces for the f-p model agree with 
the RLAPW ones. By further analyzing the f-p model, 
we have concluded that B® among the CEF parameters 
plays an important role to determined electronic proper- 
ties, consistent with the experimental results oberved in 
Celrlns and CeCoIns. 

Here it is emphasized that the system can be sensitive 
to such small energy quantities as CEF parameters, since 
the energy scale of the electron system is much reduced 
due to the large renormalization effect in heavy-fermion 



1997) 



Relativistic Band-Structure Calculations for CeTIn5(T=Ir and Co) 



11 



compound. Namely, our present conclusion is based on 
the pre- formation of heavy quasi-particles. To compelete 
the discussion, we need to treat both the effect of CEF 
and electron correlation at the same footing on the basis 
of the microscopic model, for instance, the f-p model 
with Coulomb interactions. This point will be discussed 
elsewhere in future. 

In order to make further step to quantitative discus- 
sion on supercondiuctivity, it is also necessary to ana- 
lyze the model, in which short-range Coulomb interac- 
tion terms are added to the present tight-binding model 
or more simplified /-electron model,EZP leading to a mi- 
croscopic Hamiltonian to discuss magnetism and super- 
conductivity of /-electron systems.E!r We believe that in 
this paper we show one route to arrive at the elucidation 
of the mechanism of unconventional superconductivity in 
heavy fermion materials. 
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